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Abstract 

We investigate Superconductor/Ferromagnet (S/F) hybrid structures in the dirty limit, described 
by the Usadel equations. More precisely, the oscillations of the critical temperature and critical 
current with the thickness of the ferromagnetic layers are studied. We show that spin-flip and 
spin-orbit scattering lead to the decrease of the decay length and the increase of the oscillations 
period. The critical current decay is more sensitive to these pair-breaking mechanisms than that 
of the critical temperature. These two scattering mechanisms should be taken into account to get 
a better agreement between experimental results and theoretical description. We also study the 
influence of the interface transparency on the properties of S/F structures. 

PACS numbers: 74.50. +r, 74.80.Dm, 75.30.Et 
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I. INTRODUCTION 



It has been known for quite a long time that superconductivity (S) and ferromagnetism 
(F) are two antagonistic orderings and that they can hardly coexist in a same compound. 
Although they tend to avoid each other, their mutual interaction may be studied when 
they are spatially separated, which is realized in Superconductor/Ferromagnet (S/F) hybrid 
structures (as a review, see Indeed, in such systems, superconductivity and 

ferromagnetism can influence each other through the so called 'proximity effect'. The main 
peculiarity of the proximity effect in S/F structures is the damping oscillatory behavior of 
the superconducting correlations in the F layers (while they monotonically decay in normal 
layers of Superconductor/Normal metal structures). In the dirty limit and for large exchange 
field, the characteristic lengths of the decay and oscillations are the same. They are given by 



£f = WDf/h, where Df is the diffusion coefficient in the ferromagnet and h is the exchange 
field acting on the electrons spins. This unusual proximity effect leads to several striking 
phenomena, such as the non-monotonic dependence of the critical temperature and current 
of S/F multilayers on the F layers thickness and the realization of the so-called it junction 



in S/F/S trilayers 



m 

Although the existing theory provides rather good qualitative description of the observed 
effects, there is still no complete quantitative agreement with experiments. This indicates 
that besides the exchange field, some additional pair-breaking mechanisms are present in 
the F layers. Indeed, spin-flip process is inherent to the ferromagnetic layers (because 
of magnetic impurities, spin-wave or non stoichiometric lattices) and may have dramatic 
consequences on superconductivity (contrary to non magnetic impurities that have very little 
impact). Such a pair-breaking mechanism also arises in usually used weak ferromagnetic 
alloys, because they are close to ferromagnetism disappearance and then quite favorable to 
large magnetic disorder. This can be inferred for instance from the very strong decrease of 
the critical current of S/F/S junctions as a function of the thickness of the ferromagnetic 

inn 

layer in experimental studies |4|, |5J, |f| . In such experiments, the ferromagnetic alloys used 
were Cu x Ni 1 _„ a; with x ~ 0.5, limit range of concentration for ferromagnetic properties. In 
addition, the pair destruction due to spin-orbit interaction must be taken into account as 
well. 

Though the spin-flip and spin-orbit interactions in the ferromagnetic material were in- 
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troduced in several recent theoretical studies on proximity effect in S/F hybrids, only few 
simple limiting cases were considered (we'll give appropriate references below). Thus the 
problem of quantitative description of these effects in S/F multilayered systems is still un- 
solved. In this paper we present the results of detailed theoretical study of the influence of 
spin-flip and spin-orbit scattering mechanisms on the critical temperature (T*) and critical 
current (I c ) of S/F multilayered systems. We obtain analytical and numerical solutions of 
the problem which provide the basis not only for qualitative understanding of experimen- 
tal results but also to fit the data quantitatively. In particular, we show that spin-orbit 
and spin-flip scattering mechanisms influence differently the properties of S/F structures: 
spin-orbit mechanism can destroy the T* and I c oscillations while spin-flip scattering can 
only modify them. We also report on a striking non-intuitive behavior that the critical 
temperature and current can acquire with the variation of the S/F interface transparency. 



II. LINEARIZED USADEL EQUATIONS 

A very convenient set of equations for an inhomogeneous superconductor was elaborated 
by Eilenberger j?| . However, the Eilenberger equations can be replaced by the much simpler 
Usadel equations [^] when the electron scattering free path in S/F systems is shorter than 
the superconducting length, which is often the case. These equations are nonlinear but 
can be simplified when the temperature is close to the critical temperature T c or at any 
temperature in the F layer when the transparency is low. We consider a S/F multilayered 
system where all physical quantities depend only on the coordinate x perpendicular to the 
layers. The natural choice of the spin-quantization axis is along the direction of the exchange 
field. In the general case, magnetic and spin-orbit scatterings mix up the up and down spin 
states. Therefore, two anomalous Green functions, namely F + ~ (V't'^l) anc ^ F- ~ (V'J.V't) 
are needed to describe this situation. The linearized Usadel equations may be written for 



u > as til 



in the S layers and 



w "T^ )F±s( ^ x) = A(x) ' (1) 



D f d 2 , ., 1 2 \ „ , x f 1 1 



w - -f-^ ± ih + - + - ) F±f (w, x) + I — - - ) (F ±/ ( U ,x) - F Tf ( w , x)) = 0, (2) 
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in the F layers, where A (x) is the superconducting order parameter, h is the exchange field, 
D s (Df) is the diffusion coefficient in the S(F) layer and uj are the Matsubara frequencies, 
uj = 27rTm-|-~). The parameter t so is the spin-orbit scattering time. The magnetic 
scattering times are r z = T2S 2 / (S 2 ) and t x = r 2 S' 2 / (S%). The rate r^ 1 is proportional to 
the square of the exchange interaction potential (the notations are the same as in [10]). Note 



that the microscopical Green functions theory of superconductors wit 



1 magnetic impurities 



and spin-orbit scattering was proposed by Abrikosov and Gor'kov [11 1. 

Besides, the usually used ferromagnets present a strong uniaxial anisotropy. In that 
case, the perpendicular fluctuations of the exchange field are suppressed, that is r" 1 ~ 0. 
Therefore, henceforth, r z will be noted as the magnetic scattering time r m . The spin-flip 
scattering is now simply incorporated replacing uj by uj + — in the standard Usadel equation, 
and T m is the magnetic scattering time (see for example 12j). Note that the Usadel equations 
in the F layers are not coupled anymore when = 0, and only one equation is needed. 

III. THEORETICAL DESCRIPTION OF THE NON-MONOTONIC DEPEN- 
DENCE OF T*. 

The common feature of all the S/F bilayered and multilayered heterostructures is the non 
monotonic evolution of the critical temperature T* with the thickness of the ferromagnetic 

n 

ayer. This behavior was first predicted by Buzdin and Kuprianov [13] and Radovic et al 



14j, and was since then intensively studied both theoretically and experimentally (as a 
review, see [l|). The presence of magnetic scattering can result in an additional decrease 
of the transition temperature (on the contrary, nonmagnetic impurities do not affect the 
transition temperature). Note that the theoretical description of T* for ferroma gnet ic layers 
with spin-orbit scattering was performed by Demler et al. and Oh et al. jiol] . In the 
present section, we therefore mainly focus on the influence of spin-flip process and neglect 
spin-orbit scattering. We report here on the influence of the spin-flip scattering on the non 
monotonic dependence of T* when the thickness of the superconducting layer is supposed to 
be small, i. e. d s -C £ s . In that case, an analytical solution may be obtained. It should also 
be underlined that the question of the spin-flip role was first addressed by Tagirov in 
the discussion of T* of SF multilayers and by Fal'ko et al. Q] in the study of the resistance 
of a diffusive F/N junction. 
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FIG. 1: Geometry of the considered system. 



A. Influence of spin-flip scattering on T* 



We consider a S/F multilayered system with F layers of thickness df and parallel mag- 
netization directions and S layers of thickness d s , see Fig. 1 (this system is also equivalent 
to a S/F bilayer of thicknesses df/2 and d s /2 respectively). Further, we assume that the SF 
interfaces are not 'magnetically-active', i.e. there is no rotation of the quasipartice spin at 



the interfaces as considered in Refs 



• Il9, 



20 



21 
I 



22j. Under these conditions, long-range spin- 



triplet superconductivity does not appear [3|]. Discussion of the role of spin-flip scattering 
in a ferromagnet in combination with magnetically- active SF interfaces and/or noncollinear 
magnetizations requires separate study. 

The anomalous Green's function F s varies a little in the S layer and may be approximated 



by a simple expansion up to the second order x 2 (see for example 

F s (x,uj) = F Q 



2311 and 



m) 



1 fiu 2 
1 X 

2 



(3) 



where F is the value of the anomalous Green's function at the center of the S layer. More- 
over, in that case, the spatial variation of the pair potential A(x) can be neglected A(x) ~ A. 
It follows from Eq. that Fq = where r" 1 = Rf(3 u is the complex pair-breaking 

parameter. 

The Usadel equation in the S layer is completed by the self-consistency equation 



A(x) In 



/Afar) 

it i \ I 



F± a (x,u) =0 



(4) 



where T c is the bare transition temperature of the superconducting layer in the absence of 
the proximity effect. Hence, this self-consistent equation gives for T* 

1 



T 1 * 

ln^ 



(5) 



If the temperature variation is small (r s 1 <C T c ), (JSJ) becomes 



Zk_^ = JL Re ( -i). (6 ) 



The boundary conditions for the linearized Usadel equation are 

'dFA a n fdF f 



^ X ' ±d s /2 G s \0x y ±ds j 2 



F a (dj2) = F f (d a /2)-£ nlB (?p-) , (7) 

\° X J d s /2 

F s (-4/2) = F f (-4/2) + i nlB (8) 



with a n (a s ) the conductivity of the F(S) layer, £ n = J -^v an< ^ ^ B = "^f" 2 " * s ^ e interface 



transparency, related to the S/F resistance per unit area -Rft.They lead to the determination 
of the pair-breaking parameter t~ x 

_! Dscxn F' f (d s /2)/F f (d s /2) 
Ts d s a s l- i nlB F' } {d s /2)/F f {d s /2) ' [ } 

Next, the resolution of the Usadel equation in the F layers (J2J) with symmetry consider- 
ation gives rise to the expression of the anomalous Green's function in the F layer 

F f (x, u > 0) = A cosh [k (x - 4/2 - d f /2)\ , (10) 

in the phase and 

F f {x, to > 0) = B sinh [k {x - 4/2 - d f /2)\ , (11) 

in the 7r phase. If T c < r" 1 , h, we may neglect the Matsubara frequencies in the expression 
for k that becomes 

V2 n 1 ■ 1 

k = —Vi + a = — + i—, (12) 

4/ 4/i 4/2 

with £f — J^jj- and a = The two parameters £ji and £^ are respectively the decay 
characteristic length and the oscillations period, and may be written as 

Oi = r== . (13) 

1 + a 2 + a 



« 2 — a 



For t^ 1 = 0, Eqs. (|T3j) and ((Hj) reduce to £/i = £/2 = £/■ As expected, it is found that the 
decay length and oscillations period are the same in absence of spin-flip scattering. 
The pair breaking parameter may be determined and it reads 

gtanh (qdf/2) 

r^(u > 0) = To" 1 V J , , (15) 

1 + 7gtanh ( qdf/2 

in the phase and 



gcoth ( qdf/2 

r-i(u > 0) = Tq- 1 V / (16) 

1 + 7gcoth iqdf/2) 

in the ir phase, where r _1 = ^^7^, 7 = jjIb, q = Hf and d f = 

First, if the interface is supposed to be transparent, analytical expressions of the variation 
of the temperature may be found if C T c <C 1: 

4r , t t ^ = 1 asinh(o^) - 6sin(ad^) ^ 

ft 2 cosh 2 (a<ij) cos 2 (bdf) + sin 2 (bdf) sinh 2 (a<i/) 

(18) 

4tq^ t ^ = 1 ^asinhfaj/) + frsin(kT / ) ^ 

ft 2 cosh 2 (ac?/) sin 2 (6<ij) + cos 2 (bdf) sinh 2 (a<i/) 



where two dimensionless parameters have been introduced, namely a — & — \/ \J\ + 0? + a 
and b — 4*- — V \/l + a 2 — a, so that q = a + ib. 

In the general case however a numerical analysis has to be performed. 

The ratio of the characteristic lengths 

C/2 a/ Vl + « 2 + a 

clearly shows that the magnetic scattering decreases the decay length and increases the 
oscillation period. If r m x >> h, £/i can become much smaller than £y 2 - Besides, the 
decrease of £f% makes the observation of the oscillations more difficult. 

The evolution of the critical temperature without and with spin-flip scattering is given in 
Fig. 2. The phase that really occurs is the one with higher critical temperature. It is seen 
that £fi decreases in presence of spin-flip while £/2 increases. 
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FIG. 2: Influence of the spin-nip scattering on the evolution of critical temperature as a function 
of the ferromagnetic layer thickness. 

B. Influence of the interface transparency parameter on T* 

The influence of the interface transparenc y p arameter 7^ on T* of SF bilayers was studied 
before by Aarts et al. Fominov et al. |22[ and Tagirov |Tj|. Here we extend this dis- 
cussion taking into account spin-flip scattering in the F-layer and considering broader range 
of interface transparencies and demonstrate that new effects take place in this situation. 

The intriguing evolution of T*° with the interface transparency parameter 7 must be 
underlined (see Fig. 3). If a = 0, there is no magnetic scattering and it could intuitively be 
believed that, the higher the barrier, the less is the influence of the proximity effect on the 
S layer and therefore, T* (7 ^> 1) > T* (7 ~ 1). However, it can be seen from Fig. 3 that 
the critical temperature is a decreasing function of the interface transparency parameter 7 
for a small thickness of the F layer. This counter-intuitive behavior can be qualitatively 
understood for a S/F bilayer. The probability for a Cooper pair to leave the S layer is 
smaller for a low transparent interface (7 ^> 1). Nevertheless, the probability for this pair 
to come back again in the S layer is much higher for a transparent interface. Indeed, when 
the F layer is thin enough, the reflection of the Cooper pair at the other interface of the 
F layer allows the pair to cross again the first interface, which is easier when 7 is small. 
Consequently, the staying time in the F layer increases with the increase of the barrier, 
and when this time becomes higher than the coherence time of the Cooper pair, the pair 
is destroyed, leading to a weakened superconductivity. Therefore, the critical temperature 
decreases with the barrier in that case. On the other hand, if df increases, the Cooper pair 
is hardly reflected by the external interface of the F layer whatever the value of 7 is and the 
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FIG. 3: Evolution of the critical temperature versus the interface transparency parameter for a 
S/F bilayer. 

critical temperature is expected to increase with the barrier. 

Besides, an analytical approach using formula © leads to the determination of the critical 
temperature when df < 1. In that case, the pair breaking parameter in the phase becomes 



Re (r- 1 ) 



7 + if 



(j + If 



+ 



2 ' 



(21) 



Two behaviors emerge whether jdf » 1 or not. Indeed, if 7 > df , the critical tempera- 
ture evolution is described by 



C 



7T T, 



-1 



■ydf>l 



4T C 7 



(22) 



which is an increasing function of 7. On the contrary, if 7 < df , the critical temperature 
decreases when 7 increases, following 



— =1 d f r n 7. 

T I ~ 4T 

- tc /7d / <l 



(23) 



The critical temperature evolution with 7 presented in Fig. 3 is therefore understood. 

In presence of small magnetic scattering, the critical temperature evolution remains qual- 
itatively the same. As shown in Fig. 3, it appears that superconductivity may be less de- 
stroyed when there is spin-flip scattering. In that case, the pair breaking parameter becomes 



Re (tJ) 



7 + 



+ ^ 



(24) 
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The critical temperature is described by 



+ ( 25 ) 



when dfj > 1, and magnetic scattering leads to a slight enhancement of the transition 
temperature. On the contrary, ii^df < 1, 

'£)_ = (£<^)_ -f^-y (2 6 ) 

and spin-flip implies the decrease of the transition temperature. 



C. Influence of spin-orbit scattering and 'perpendicular' spin-flip 

Let us now consider briefly the general case, with spin-orbit and/or perpendicular spin-flip 
scattering. An additional parameter has to be introduced, namely 

°x=i(l-^), (27) 

and the parameter a now becomes 



In that case, expressions (J15|) and (|16|) are modified. In the phase, gtanh [qdf J is replaced 
by 

g*tanh (q*df/2^ — gtanh (qdf/2 



gtanh (qd f /2 +(3 V TT^J^V S ( 29 ) 

P + 



l+7g*tanh(q*(i / /2) 
1+79 tanh (qdf /2 

where q becomes 



q 2 = 2[^ + a-a L +isJl-a\) , (30) 



a 2 * 



and /3 is 

n. /'v'l 

= -a± = • (31) 

1 + a/1 - aj_ 

In the 7T phase, gcoth (qdf/2) is replaced by 



<fcoth (q*d f /2^J - gcoth (V//2 



,00^(^/2) + p A ^gg^ (32) 
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FIG. 4: Monotonic evolution of the critical temperature for an transparent interface, and for 
ro = 22 (cgs), a = 1 and a± = —3/2. 

Therefore, the influence of 'perpendicular' spin-flip scattering and spin-orbit scattering is 
quite similar to the influence of 'parallel' spin-flip processes, in the sense that it also implies 
the decrease of the decaying length and the increase of the oscillations period. However, a 
special situation arises when a± > 1. Then, the oscillations of the Cooper pair wave function 
are completely destroyed. Similar conclusion for spin-orbit mechanism was obtained in ^| . 
In fact, the influence of the 'perpendicular' magnetic scattering is analogous to the spin-orbit 
scattering (see Fig. 4). Probably the role of 'perpendicular' spin-flip or spin-orbit scattering 
is important for the understanding of experimental results where no oscillation of the critical 
temperature was detected. Besides, note that the critical temperature oscillations can not 
disappear when there is only 'parallel' spin-flip. 



IV. BEHAVIOR OF THE CRITICAL CURRENT I c 

A Josephson junction may be realized with S/F/S sandwiches in which the weak link 
between the two superconductors is ensured by the ferromagnetic layer. The supercurrent 
J s (ip) flowing across the structure can be expressed as J s ((f) = I c sin (<p), where I c is the 
critical current and ip stands for the phase difference between the two superconducting layers. 
A standard junction has at equilibrium J c > and ip = 0, and therefore, no current exists. 
It may appear however that I c becomes negative, which implies that the equilibrium phase 
difference is tp = ir and the ground state undergoes a 7r phase shift. This so called V junction' 
was first predicted for S/F/S structures in the clean limit [28.], and later in the more realistic 
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case of the diffusive limit [29f . The critical current I c in S/F/S junctions is controlled by the 
exchange field in the ferromagnet, the interface transparency parameter 75, the spin-flip and 
spin-orbit scattering rates. In this section we will discuss the influence of these parameters 
on I c , paying the main attention to the role of spin-orbit and magnetic scatterings. Besides, 
note that the spin-orbit effect on the critical current has already been studied for S/F 
bilayers separated by an insulator with thin F layers 3(| and for the FSFSF structure 
with noncollinear magnetizations of F- layers We will consider a symmetric S/F/S 

structure with the F layer having thickness df and assume that the dirty limit conditions 
are fulfilled for S and F materials. We will assume that the rigid boundary condition 
7b >> m.m{a n ^s/o'sCn, 1} is fulfilled, when the suppression of superconductivity in S can 
be neglected and the linearized Usadel equations © can be used. The case of transparent 
interfaces, 7# = 0, will be considered at the end. 



A. General expression for I c 

To derive the general expression for I c we should supplement (j2J) by the boundary condi- 
tions (25! at S/F interfaces (x = ±d/) 



lB^F ±f (±d f ) = ±G S 



A(±d f ) 



10 



F±f(±d f ) 



sign{uj), 



(33) 



00 



-., A(±d f ) = |A| exp {±iy>/2} . 



>lu 2 + \A\ z 

The solution of the boundary problem ©, (J3~3*|) has the form 

F + f = ^((l+r^) {A u+ cosh [k + x] + B u + sinh [k + x]}+(l-r)u) {Aj- cosh [k-x] + sinh [k_x]}), 

(34) 



F-f = -((1—riu) {A u+ cosh [k + x] + B^+ sinh [k + x\}-(l+r] u ) {Aj- cosh [k-x] + sinh [k^x]}), 



where the coefficients A w ± and are given B^ in the Appendix, while 



+ v /q4 3 T , V2 fp— - 7= - 

, / n — 7 k± = —y + a so + a± V a z so - 1. 



ihsign((jj) 



it 



(35) 



(36) 
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The parameter a so is defined by a so = l/(hr so ), a = l/(hr m ) and u = uo/h. Note that in the 
general case 'perpendicular' spin-flip has to be added in a so and a. The general expression 
for the supercurrent is 



• rri °U 

ml a n 

J s - — — 2^ 



4e 



W = — 0O,<7=± 



dx dx 



(37) 



where F±f(x,u) = F±* (x, — uj) . Substituting (|3*4*j) and (|3*5|) into the above expression and 
taking into account the symmetry relations (given in the Appendix), we get Js = I c smip. 
For large spin-orbit scattering, i. e. t~J- > h, the critical current I c equals to 



27rT c a n T c ^ 2 (l + r^) 1 + + 



where 



inUl k ± + G 2 s) sinh [M/] + 2G slB £ n k ± cosh [^d/] ' 
At larger exchange field h > r," 1 , the critical current becomes 



eL 



T 



27iT c a n £ f T c 



oj=0 



Re 



(a + z6) 


f~2 

V OL so — 1 




((7bM 2 + G|) sinh 


fed/] + 2Gs^Bink cosh [fcdy] 



where k — q/^f and q = a + ib, with 

a = y|£5| + a so + a + + a ) 2 + 2(p| + a)a so + 1, 



(38) 



(39) 



(40) 



(41) 



1 - a 2 



(42) 



\lu\ + a so + a + a/(M + a) 2 + 2(jtJj + a)a so + 1 
I n ap pro p riate limits, expressions (|38j h (|40j ) transform into the previously obtained results 



in appror 

Q, Q, 



24j,j32j. This approach is valid if F w ± C 1 or 

a A 1 

7b > 7, 



+ 5 2 7rT c min{l,d//£ n }' 



(43) 



As follows from (|4()jl - (|42jl . in the case of strong spin-orbit scattering r^ 1 > h, the critical 
current decays monotonically with the increase of df with two decay lengths k± x defined by 
Eq.(|36J). It is seen from (jHUjl that in the limit h — > 0, the parameter rj w — > oo. As a result, the 
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contribution to the critical current in ()38j) comes only from k_ component with the length 
scale kZ 1 which describes the case of an S/N/S junction in which spin-orbit scattering does 
not influence I c . With the increase of h, the contribution to J c from the faster decaying 
k + component (k + > kJ) also increases and the difference between k + and k_ decreases. 
Finally, when h = r" 1 , both scales coincide, k + = k_, and the components k + ,k_ provide 
equal contributions to I c . 

For relatively weak spin-orbit scattering, r" 1 < h, the dependence I c (df) follows the 
damped oscillation law, when two length scales £/2 can be introduced describing respec- 
tively the decay and the oscillation period of I c (df). In this section, we will concentrate on 
the case t~J- < h. The scales £/ 2 ar e related to a, b and will be discussed in detail below 
in different limits. 



1. Limit of small F layer thickness and large jb- 



In the limit of small thickness df <C and large interface transparency parameter 
7_b 3> df/£ n , { r KT c Gsin)/{df{uj + TZ n l + t~q) we can neglect the terms G 2 S and GsjB^nkikdf) 2 
in the denominator of Eq. (}4T)j) and with the accuracy of better than ((i//£/i) 3 get 



AixT r . 



jdfT t 



oj=0 



a 2 + v 2 + + v 2 ) 



where Rn = 2Rb is the normal junction resistance, 7 = jB^n/Cf, uj = u/h and 



If additionally 



= lu + a so + a, 



IB > 



fix = tt + 2Gc 



j B hdf' 



then 



d f 3 [(nT c /h + a so + a) 2 + 1 - a 2 J] h 

rp OO 

-T 



£/ J T ^ A 2 + 2« so d/ A A 

o — ~~~ t cLllll — — 

6^ 2T C 2T 



eRiyI c 

4vrT c tydf \ T c ^ u 2 + A 2 (ft 2 + u 



(44) 



(45) 



For h ^> 7rT c , Tm 1 , r^ 1 and T <T C , the sum in (J43|) can be calculated by transforming from 
summation into integration over u resulting in 



eR N I c A £f 



AixT c 2T c ^d f 



2A h 
— m — 

llTT A 



hr s , 



+ 



dj_ 



(46) 
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From ()46|) it follows that the transformation from to it junction may occur when the 
thickness df exceeds some critical value d n 



For t^" 1 ,^ 1 <C h this result agrees with |32j]. It is interesting to note that in this range of 
parameters, the condition (}4Tj) of the transition to 7r state does not depend on the properties 
of the interfaces. As is also seen from Eq. ()47j). both spin-orbit and spin-flip scattering 
increase the thickness d n corresponding to the first to it transition. 

If'jdf/C.f < 1, then the term GsCf /{idf) m is not small. As a result, the transition 
to a 7r state should depend on the properties of interfaces and occurs at h larger than the 
critical value h n oc 7rT c (£ ra /<ij) following from Eq. (}4l)j) . 

Therefore, in the limit of small df <C £/i and large 7# ^> 1 the transition from to it 
junction exists only if the exchange energy h > h w oc 7iT c (^ n /df) sufficiently exceeds ttT c . 
The smaller the interface transparency parameter js, the larger should be the exchange 
energy h. The possibility of the 0-7r transition in junctions with small thickness df in the 
case of large 7b is related to the multiple scattering at the boundaries. As a result, the 
electrons spend more time in the F layer and its influence on superconductivity is enhanced. 
This is the manifestation of the same mechanism that leads in S/F bilayers to the critical 
temperature decrease with the increase of the interface transparency parameter (see section 




(47) 



II B). 



2. Limit of large F layer thickness. 



In the limit of large df S> £/i for the critical current from ()4()|) - (021), we have 



el c 



2T 



oo 



A 2 G 2 Q usin(bdf/C,f) + v cos(bdf/^f) 
((G 5 + a7) 2 + 7 2 6 2 ) 2 



{-^}, (48, 



2nT c a. 



E 



exp 



n 



where the coefficients u and v are defined by 





G 2 S (a + 7(a 2 + 6 2 ) (2G S + cq)) , 



2 



(49) 
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If additionally min {h, t s J-, r m 1 } 3> nT c , then both a and 6 may be considered as inde- 
pendent on Matsubara frequencies, since the sum in (|4*Ujl converges at uj « 7rT c . In this case 



27rT c o- n 

where 



^ ■ df ^ df 
hi sin — — + L 2 cos — — 

f/2 f/2 



exp^-^|>, (50) 



_ 2T ^ A 2 G| 
1 T c ^ o ^ ((Gs + ^ )2 + &2 ^ 2) 2' 



2T ^ A 2 G| 

The two characteristic length scales are given by the following expressions 



Oi=W— 7 7 2^0 T f ' (53) 

a + a so + v a + 2aa so + 1 



, /a + a so + \/a 2 + 2aa so + 1 
2 = W f T , (54) 



so 



which generalize Eqs. (fTTjj) and ()14|) for the case of the presence of the spin-orbit scattering 
in a ferromagnet. One can see that with the increase of both scattering rates r" 1 and r^ 1 
the decay length £ fx decreases, while the oscillation period £ f 2 increases. 

For a weak exchange field h << ixT c and sufficiently high temperatures t\T >> 
{t~J-, r" 1 , h}, with h > r" 1 , only the first term with n = in (j4T)j) is important and we 
have 

eC f I c _2A 2 G Z lS m(d f /£ f2 ) + i: 2 cos(d f /£ f2 ) \ d f 



2ixT c a n tt 2 TT c (( Go + a^f + 6^) 2 6XP \ £/i / 



where Gq — kT/ \J (vtT) 2 + A 2 and Ei and E 2 become 



2T A 2 Gp 

T c (vrT)2 ((Go + a ^ )2 + 6 2~ 2 



Sl " -^TZ^— . ° ,2-2,2' (56) 



2T A 2 G 2 



v 



T c (7rT) 2 (( G + a0 7)2 + 6^2)2- ( 57 ) 

The parameters a , &o, «0j are obtained replacing uj by tcT in expressions (jUJ, (f4*2|) . (|4T)|) . 
The two characteristic lengths may be written in that case as 



Zn = Zf\ / > (58) 

7rT/ /i + a + a so + a/ (vrT/ ft, + a) 2 + 2{ixT + a)a so + 1 
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'nT/h + a + a so + y/(nT/h + a) 2 + 2(ttT + a)a so + 1 

l-ai 



& = Sf\l 1 — 3 • ( 59 ) 



so 



For TgJ-jT^ = 0, Eqs. (|H%|) . (jH^j) reduce to a simple expression £* 12 = ^ n 1 [(/i 2 /(7rT c ) 2 + 
(T/Tc) 2 ) 1 / 2 ± (T/Tc)] 1 / 2 |4| which describes the temperature variations of both length scales. 



More precisely, in the considered limits, ^T 1 ~ £ s ~ 1 yj2T/T c , while ~ £7 1 1 /i/(27Ta/TT c ) << 
£7 X . One can see that the scattering rates t^ 1 ,?"" 1 make the T variation of £/i t 2 weaker. 
From 15()n . ()55jl . one can derive the condition for J c = and obtain 

df n = df, = 7rn — arctan(— -), n = 0,1,2,... (60) 

when the transitions between 0— and 7r-states occur. The position of the first zero, dfi, 
depends both on the material parameters of the ferromagnetic layer and the properties 
of the interfaces and superconducting electrodes, while the distance between the zeros is 
the function only of and therefore depends only on the transport parameters of the 
ferromagnetic material. 

From the structure of coefficients u and v (see Eq. (|4*9*|) and (}60|) ). it follows that in the 
limit of small 7# (F±f = ^^_ a2 exp (±iip/2) at S/F interfaces) and min {h, t~J-, t" 1 } 3> 

^ = irk - arctan( 2& ~ av/ < ^), fc = 1, 2, ... (61) 

In particular, it follows from ()61|) that for small spin-orbit and spin-flip scatterings h ^> 



r- 1 r- 1 

'so I'm i 



d< 



71 



-g™*k-^ (62) 

and the well known result df% m (37r/4)£/2 for the first critical thickness for to tt state is 
reproduced. This df\ value approximately satisfies the condition of validity of the large df 
approximation considered in this section. 

It follows from (JBTj) that the critical thickness df\ increases with r m x and (see also 
numerical results below). 

An increase of 7^ results in the suppression of the magnitude of F±f(±df) ~ A^I + jb) -1 
near the S/F interfaces, which leads to the decrease of df\. Formally it is due to the fact 
that with an increase of 75, the coefficient b in (149(1 and, hence, Ei, may change its sign 
resulting in the existence of the solution df of (|5U|) for n — 0. This solution corresponds 
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to the thickness range in which simple large df approximation (J4*%j) is no valid anymore. 
This fact is in the full agreement with our consideration performed for the limit of small df. 
Namely, it follows from (|47|) that at large 7^ and h we have 



df! = v^/y ^ ln ^ + a + 3 «so < £fi- (63) 

Note that in actual experimental situations, when the approximation df ^> £/i is not 
fulfilled, simple expressions (JoTf)) . (|53|) . (J52J), (fo"%|). (jo^|) . fjo^|) are not valid and to analyze the 
data it is convenient to introduce the effective decay length tffi 

#/ = _4_]n^, (64) 

"ml w m 2 lcm2 

where d m \^ are the thicknesses at which the first (I cm 2) an d second {I ana) maxima of I c (df) 
occur. 

Below we will focus on the influence of jb, t~J-, and r" 1 on the critical current. 
B. Influence of interface transparency parameter j B on I c 

Consider first the simplest case of vanishing r,," 1 , r" 1 . In this limit, we immediately 
deduce from lj5njl - (|32jl that 



27rT c (7 n T c ^ ^2 5 Re ( (( 7B ^ n A;)2 + G|) S i n h [fcd,] + 2Ggy B i n k cosh [fed/] ' ' (65) 



aj=0 

where the coefficients a and b become 



Z\ + J\u\ 2 + l, (66) 



£| + yp| 2 + l. (67) 

The dependence of the critical current as a function of the F layer thickness calculated from 
(|65|) for T = 0.5T C , and /i = 3tcT c for different values of the interface transparency parameter 
7s is presented in Fig. 5. Note that for small 7# Eq . (J65j) is formally not applicable, since is 
was derived under the condition (|4"3*|) of sufficiently large 75 . Therefore, for 7^ =0, -ys =1 
and 7b =5 we have added for comparison the corresponding curves calculated numerically 
by direct solution of the Usadel equations for arbitrary j B , which show that Eq. (j63j) provides 
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FIG. 5: (Color on line) Influence of the interface transparency parameter 7^ on the thickness 
dependence of the critical current in S/F/S junction for a = a so = 0, h = 3irT c and ^- = 0.5. 
Open circles: the results of exact numerical calculations for jb =0, jb =1 and 75 =5 (from top 
to bottom). 

reasonable approximation even in the small jb range. It is clearly seen from Fig. 5 that 
with an increase of jb, the position of the first zero dfi is shifted into the region df < 
As discussed above, the increase of ■jb results in the suppression of F±f(±df) in the F layer, 
such that ReF±f(x) changes sign in the F layer center at smaller df . 

The results plotted in Fig. 5 make it possible to estimate the upper limit for the I c Rn 
product of the SFS junction in a 7r-state. This upper limit can be achieved in the case of 
highly transparent interfaces (7_b=0) and in the absence of spin-flip and spin-orbit scattering. 
As follows from Fig. 5, IcR™™ ~ 0.l7rT c that provides I c R]} ax « 250 /iV in the case of Nb 
electrodes (T C =9K). 

At large F layer thickness, I c (df) is determined by Eqs.(|5D|). (j53|) . If spin-flip and spin- 
orbit scattering are negligible, the ratio £/i/£/2 = h/(T + a/ (tiT) 2 + h 2 ) follows from (p)%j) . 
(|59|) and depends only on h and T. For typical ferro magnets h > 7rT c and £/i approximately 
equals to £/2- However, if the spin-flip scattering and spin-orbit rates become relatively large 
T m 1 5 r so 1 ~ the situation may change drastically. Consider first the influence of spin- flip 
scattering on I c . 
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FIG. 6: (Color on line) Influence of the spin-flip scattering parameter a on the thickness dependence 
of the critical current in S/F/S junction for a so = 0, h = 3irT c , 75 = 10, and = 0.5. 

C. Influence of spin-flip scattering on I c 

If spin-flip scattering is not negligible, then the ratio of the characteristic lengths in the 
decaying solution (J5UJ) . (|55|) becomes 

^ = r ( 68 ) 

and for strong spin- flip scattering r m x > /i the decay length may become much smaller 
than the oscillation period £/ 2 - This results in the much stronger decrease of I c versus df in 
S/F/S junctions. 

The evolution of I c for different values of r" 1 calculated from (|40|) for T = 0.5T C , /i = 3nT c 
and 7b = 10 is given in Fig. 6. 

One can see that with increasing a, the critical thickness dfi of the first — n crossover 
shifts to larger values of df. 

D. Influence of spin-orbit scattering on I c 

If spin-orbit scattering is not negligible, then the ratio of the characteristic lengths in the 
decaying solution (jSOJ), becomes 



^ = — V * ~ = — (69) 

O2 + a so + a/w 2 + 2 |w| a so + 1) 
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FIG. 7: (Color on line) Influence of the spin-orbit scattering on the thickness dependence of the 
critical current in S/F/S junction for a = 0, h = 3ttT c , 75 = 10, and = 0.5. 

Eq.(|69j) shows that the difference between the decaying £/i and the oscillating length £p in- 
creases with t~ 1 even faster than for the case of spin-flip scattering. Moreover, the transition 
to monotonically decaying solution takes place at r" 1 — > h and £/ 2 — 00. 

The dependence of critical current on df calculated from f|4*0|) for T = 0.5T C , h = 3ttT c 
and 7s = 10 is shown in Fig. 7. One can see that the oscillation period increases strongly 
with increasing a so and diverges when a so = 1. At the same time, dfi shifts to larger values 
of d f . 

V. CRITICAL CURRENT OF LONG JUNCTIONS WITH TRANSPARENT IN- 
TERFACES AT ARBITRARY TEMPERATURES. 

Let us now consider a S/F bilayer with a transparent interface. The complete nonlinear 
Usadel equation in the F layer has to be employed. For the easy axis of the ferromagnet 
and using the usual parametrization of the normal and anomalous functions Gf = cos#/ 
and Ff = sin8f, it may be written in the absence of spin-orbit scattering as 

P ° 2 n ( M + lhsi 9<^) cosflA 

t»w 6 ' ~ { — w e — + ) 8111 9f - °' (70) 

Note that Eq.([70|) transforms into (0) in the limit of small 9/ < 1. 

For S/F/S junctions, this equation may be used and an analytical solution found if 
df 3> £fi. In that case, the decay of the Cooper pairs wave function in first approximation 
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occurs independently near each interface []] . It can therefore be treated separately enough to 
consider the behavior of the anomalous Green's function near each S/F interface, assuming 
that the F layer thickness is infinite. 

For one interface (x = —df/2), a first integral of (J7U|) leads to 



U dx 2 



|a>| + ihsign(u) . Of 



1 + 



cos 



(71) 



r m (\uj\ + ihsign(uj)) 2 
where the boundary conditions Of {x — > oo) = has been used. Further integration in (j71j) 
gives 



1 -e 2 sm 2 °4] 



cos y 



g exp 



1 — e 2 sin 2 -y ) ( ' ( ^ 



-2q 



(df/2 + x) 
if 



(72) 



where k = y/2/^f\/ + isign(u) + a) = g/£/, and e 2 = a/ (\uj\ + isign{u) + a) . 

The integration constant go in (f72"j) should be determined from the boundary conditions 
at S/F interfaces. For simplicity we will assume that the rigid boundary conditions are valid 
at x = —df/2, then 



0f(—df/2) = arctan • 



(73) 



and from (J72J), (|73^1 we finally get 

£o = 



e 2 ) T (n) 



and 



^/(l-s 2 )F(n) + l + l 



(A/(2tKT)) s 



i 2 ' 



n + 1/2 + J {n + 1/2) 2 + (A/ (2ttT)) 2 



2 ■ 



where f2 = y oo 2 + |A| . Next, may be determined. 

The anomalous Green's function at the center of the F layer in the S/F/S junction may 

be taken as the superposition of the two decaying functions jsJ] taking into account the 

phase difference in each superconducting electrode, 
4 







f 



9o exp 



(df/2 + x) .<p\ \ (x-df/2) .ip 



(74) 



As a result, the current-phase relation is sinusoidal and the critical current becomes 



LR 



c-n-V 



64?r- 



djT 



( 



Re 



T (n) q 
l-e 2 )f(n) + l + l 



exp(— q 



df, 

if' 



(75) 
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The critical current is proportional to the small factor exp(—qdf/£ n ). The terms neglected 
in our approach are much smaller and are of the order of exp(— 2g<i//£ n ). Therefore, they 
give a tiny second harmonic term in the current-phase relation. 

t should be underlined that this expression coincides with the one previously obtained 
in 2^1 in the limit r m x — > (no magnetic scattering) and for T c << h. 

In addition, note that in the limit of vanishing magnetic scattering, the temperature 
dependence of the critical thickness dfk of F layer when I c = may appear only through the 
variation of k with uj (where k ~ a/2 (\uj\ + isign(uo)) /£/ in this case). As the characteristic 
range of the variation of the Matsubara frequencies in the sum in ()75j) is u ~ 7rT c , then in 
the limit h ^> 7rT c , t" 1 , the dependence of dfk on T will be weak. However, when the spin 
scattering is not weak, i.e. r m x > h, another mechanism of the temperature dependence of 
the dfk emerges due to the temperature dependent term r m 1 cos^/ in the Usadel equation 
(J7UJ), or, in another words, due to the complex uo dependent function 1 — e 2 . For a strong 
ferromagnet, h 3> 7rT c , the latter mechanism may become stronger than the 'thermal' one 
related to k(u). It is not difficult to take into account both mechanisms in the numerical 
calculation of the sum (f7H|) . 

As an illustration, we present in Fig. 8 the theoretical fit of the experimental data for 
NbCuo.52Nio.48Nb junctions by Sellier et al. [35]. making use of Eq. fj75|) valid in the limit 



of small interface resistances (75 = 0). 
while the experimental value presented in 



?he F layer thickness used for fit is df = 18 nm 
] is 19 nm. The difference may be explained by 
the uncertainty in F layers determination that is around 1 nm, and may even increase due 
to the presence of a magnetically dead layer near S/F interface. Taking this in mind, the 
theoretical description of the critical current temperature dependence can be considered as 
rather satisfactory. Besides, Houzet et al. j^J ha ve p erformed numerical calculations and got 
also good fit of another experimental curve of 35j] (df = 17 nm) for the set of parameters 
which are in the same range that we have used. 

As another application of the formalism, we present in Fig. 9 the theoretical fit of the 
experimental data of Ryazanov et al. for NbCuo.47Nio.53Nb junctions with df = 22 
nm. Good agreement is achieved assuming h = 650 K, a = 1.35 and df = 21 nm. As 
in the previous fit, small difference in df may be explained by the uncertainty in F layers 
determination. 
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FIG. 8: (Color on line) Fit to the experimental data from Ref. 35] for the critical current in a 
NbCuo.52Nio.48Nb junction. The fitting parameters are: h = 22QK and a = 3. 




FIG. 9: (Color on line) Fit to the experimental data from Ref. 36] for the critical current in a 
NbCuo.47Nio.53Nb junction. The fitting parameters are: h = 650JT and a = 1.35. 

Note that the rather complicated expression for the critical current ()75j) simplifies near 
T c and may be written as (for T c «h) 



LR 



A 2 7id f 1 . dt . ( df 

— 7TT exp ( — —r— ) sin * + ~r~ 

2eT c e /2 cos(*) Pl C/i V C/2 



(76) 



where \l/ is defined by tan(\I/) 



in 
in 



The damping oscillatory behavior of the critical 



current is therefore retrieved, and the simplicity of the previous expression makes it useful 
for theoretical description of the evolution of the critical current with the thickness of the 
F layer. 



24 



VI. CONCLUSION 



We have made detailed theoretical investigation of the influence spin-dependent scattering 
in the ferromagnet on the critical temperature of S/F multilayered systems and the critical 
current of S/F/S Josephson junctions. More precisely, we have demonstrated that spin- 
flip and spin-orbit scatterings both lead to the decrease of the decaying length and the 
increase of the oscillations period. Besides, spin-orbit scattering may be more harmful for 
superconductivity than spin-flip scattering. Indeed, the oscillations of T c and I c can be 
destroyed by spin-orbit scattering, while spin-flip scattering can only modify them. This 
allows to distinguish these pair-breaking mechanisms, which both should be taken into 
account for theoretical fits of experimental data. Moreover, the simultaneous introduction 
of the spin-flip and the spin-orbit scattering leads to important new prediction that the 
perpendicular spin flip can compensate spin-orbit and vice versa. Then the perpendicular 
magnetic scattering may restore the oscillations of T c and I c which would be otherwise absent 
if spin-orbit scattering is strong enough. 

We have also studied the influence of the interface transparency on T c and I c in S/F 
and S/F/S structures, respectively. The nonmonotonous behavior of T c with the interface 
transparency parameter 7# was predicted in the case of thin F-layers due to multiple scat- 
tering at the interfaces. We have shown that the same mechanism is responsible for the O-tc 
transition in S/F/S junctions with small thickness df in the case of large jg. For larger df 
detailed analytical and numerical study of the influence of 7b parameter is presented. 

It was predicted in Ref. 3] that in the case of noncollinear magnetic ordering the long- 
ranged triplet component of the superconducting condensate could be generated. The spin- 
orbit and the perpendicular magnetic scattering are rather harmful for this long-ranged 
triplet component. Our analysis shows that the ferromagnetic alloys like Cu^Mi-.,., where 
following our estimates the parameter 1 / hr m may exceed unity, are not suitable candidates 
for the experimental search for the triplet component. 

Finally, note that it may be interesting to study the Josephson junctions with the ferro- 
magnetic layer substituted by a paramagnetic one. Applying an external magnetic field H, it 
is possible to produce the polarization of the magnetic atoms and then generate an internal 
field h ~ xH, where \ is the paramagnetic susceptibility. The variation the external field 
allows to change the relative contribution of the scattering mechanisms and the exchange 
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field and then to modify the properties of the junction in a controllable way. In particular, 
it could provoke the transition from to tt state. 

The authors thank M. Daumens, M. Houzet, M. Kulic, C. Meyers, and V.V. Ryazanov for 
useful discussions. We also thank H. Sellier and V.V. Ryazanov for communication of their 
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VII. APPENDIX: 

The coefficients A w ± and B w ± that appear in an d may be written as 
K = 2AG s cos(y/2) ^ 

r> 2zAG 5 sin(y>/2) 



2 AG S cos(y>/2) 


M {lBink+ 


sinh [fc+rf/] + \Gs \ cosh (1 + rfc) ' 




2iAG s sm(ip/2) 




cosh [fc+c?/] + \Gs\ sinh [fc+d/]} (1 + rfc) ' 




2AGsVw cos(^/2) 


M {iBink- 


_ sinh [fc_d/] + |G S | cosh [fc_d/]} (1 + 77*) 




2iAG s VusH<P/ 2 ) 



= ^urgi^aiii^/^ , gQ . 

\u\ {^B^nk- cosh [k-df] + \G S \ sinh (1 + rfc) ' 

The symmetry relations following from (J53j) - (jSBj) are at h < 

Vu = V-„, A*_ UJ+ = A U1+ , A*_ u _ = A w - i B*_ u+ = -B u+ , B*_ u) _ = -B u _. (81) 

and if h > rZ 1 : 



(82) 



< u »7u> = -1, Pu>+ = PI- = P, 

V-u-AL u+ = Au-, A*_ L0 _ = rj w A w+1 ?7* w -£>* w+ = —Bui-, B^ w _ = —r) w B u+ . 

These relations allow us to simplify the calculations for the determination of the critical 
current. 
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